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Abstract 

The paper deals with the problem of the existence of a normal form for a nearly-integrable 
real-analytic Hamiltonian with aperiodically time-dependent perturbation decaying (slowly) 
in time. In particular, in the case of an isochronous integrable part, the system can be cast 
in an exact normal form, regardless of the properties of the frequency vector. The general 
case is treated by a suitable adaptation of the finite order normalization techniques usually 
used for Nekhoroshev arguments. The key point is that the so called “geometric part” is not 
necessary in this case. As a consequence, no hypotheses on the integrable part are required, 
apart from analyticity. 

The work, based on two different perturbative approaches developed by A.Giorgilli et ah, is a 
generalisation of the techniques used by the same authors to treat more specific aperiodically 
time-dependent problems. 
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1 Introduction 

The problem of casting an analytic nearly-integrable Hamiltonian system into normal form is 
deeply related to Poincare’s challenging probleme general de la dynamique [Poi92]. Nowadays, 
normal forms are still one of the main technical tools used to deal with the issue raised by 
Poincare in this context. 

The particular case in which the unperturbed part is supposed to be linear in the actions 
(isochronous case), already investigated by Birkhoff (and for this reason also known as the 
Birkhoff problem) [Bir27], has a peculiar interest. The first rigorous statement concerning its 
stability can be found in [Gal86]. The possibility to cast the considered Hamiltonian in normal 
form, up to some finite order^ r and to obtain, as a consequence, a stability time estimate “a la 
Nekhoroshev”, is directly related to a particularly simple small-divisors analysis: the non-resonant 

‘This research was supported by ONR Grant No. N00014-01-1-0769 and MINECO: ICMAT Severe Ochoa 
project SEV-2011-0087. 

E-mail: alessandro.fortunati@bristol.ac.uk 
E-mail: s.wiggins@bristol.ac.uk 

^It is easy to see that any attempt to consider the limit r —^ oo would imply the degeneration into a trivial 
problem, (i.e. in which the allowed perturbation size reduces to zero, see also [GG85, formula (46), Pag. 105]). 
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(Diophantine) hypothesis on the frequency vector a; of the unperturbed system is sufficient in 
order to ensure the resolvability of the (standard) homological equation arising in the normal¬ 
ization algorithm. An extensive bibliography on this problem goes beyond the purposes of this 
paper, we only mention the recent generalisations for the Planetary problem of [Pinl3] and of 
[Bam05] for infinite dimensional systems. 

ft is well known that the extension to the non-isochronous case requires a careful analysis (ge¬ 
ometric part, see [Nek77], [Nek79] and [BG86]) on the regions of the phase space in which the 
actions I are such that w = uj{I) is non-resonant (non-resonant domains). 

The problem of dealing with time-dependent perturbations without any hypothesis on the time 
dependence (e.g. periodic or quasi-periodic) has peculiar technical difficulties. After the pioneer¬ 
ing works of [Pus74] and [GZ92], the interest for this class of problems has been recently renewed 
in [Boul3], [FW14a] and subsequent papers. Basically, the novelty consists in the treatment of 
the time-dependent homological equation. A first approach consists in keeping the terms involv¬ 
ing the time derivative of the generating function (also called extra-terms) in the normal form 
and then providing a bound for them. This approach, originally suggested in [GZ92] then used 
in [FW14a], yields a normal form result for the case a of slow time dependence. This hypothesis 
provides a smallness condition for the mentioned extra-terms. Alternatively, those terms can be 
removed by including them into the homological equation, which turns out to be, in this way, 
a linear ODE in time. This has been profitably used in [FW14b], [FW15a] and in [FW15b] but 
requires (except for a particular case described in [FW15b]) an important assumption. More 
precisely, it is necessary to suppose that the perturbation, as a function of t, belongs to the class 
of summable functions over the real semi-axis^. As in (3), those functions exhibiting a (slow) 
exponential decay will be used as a paradigmatic case, ft will be shown that the consequences 
of this assumption in the isochronous case are remarkable: the normalization algorithm can be 
iterated an infinite number of times by means of a superconvergent method borrowed from KAM 
type arguments, see e.g. [Chi09]. The procedure leads to the so-called strong normal form i.e. 
in which the normalized Hamiltonian has the same form of the integrable part of the initial 
problem. Furthermore, no restrictions are imposed on w, hence flows with arbitrary frequencies 
persist in the transformed system. 

As it would be likely to expect, this phenomenon has an important consequence also in the 
non-isochronous case. The possibility to disregard the problems related to the small divisors 
implies that the well known geography of the resonances analysis, a key step of the Nekhoroshev 
theorem, is not necessary in this case and the results that can be stated are purely “analytic”, 
fn such a way, the classical assumptions on the unperturbed part of the Hamiltonian (such as 
steepness, convexity etc.), are no longer required. As a common feature with the isochronous 
case, the obtained normal form does not exhibit resonant terms, as these have been annihilated 
in the normalization by using the time-dependent homological equation. This implies that, in 
this case, the plane of fast drift (see e.g. [Gio03]) degenerates to a point. 

The paper uses in a concise but self-contained form, the tools developed in the above mentioned 
papers of the same authors, especially of [FW15b] in which the concept of “family” of canonical 
transformations parametrised by t is introduced. The proofs are entirely constructed by using 
the language and the tools of the Lie series and Lie transform methods developed by Giorgilli et 
ah, see e.g. [Gio03]. 

^We stress that this hypothesis is usually not satisfied in the case of periodic or quasi-periodic time dependence. 
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2 Setting and main results 

Consider the following nearly integrable Hamiltonian 

H{I,ip,rj,t) = h{I) + r] + ef{I,ip,t), (1) 

with G G X ¥"■ X M X M'*', where G C and e > 0 is a small parameter, which is 

the “autonomous equivalent” in the extended phase space of Hamiltonian T-L{I,ip,t) = h[I) + 

We define, for all t G := [0,+oo), the following complexified domain Pp,cr := x ¥” x Sp, 
where Gp ■= U/eG 

Ap{I) := {/ G : |/ - /| < p}, := ^ G < ^x}, 5p := {r/ G C : |9r?| < p}, 

with p, u G (0,1). For all p : ^p x T” x M"'' —>■ C, write g = 9 k{I,t)e^^'^, then define the 

Fourier norm (parametrized by t) 

llslU- E (2) 

fceZ" 

with I • Ip is the usual supremum norm over Qp and \ k\ := \ ^l\- ^ x T” x —>■ C” 

we shall set ||r(;||^^ := Yl?=i ll^dlpo- standard framework (see eg. [BGGS84]) is the space 
(tp^cr, of continuous functions on Qp x T”, holomorphic in its interior for some p,a and real on 
G X T” for alP t G We shall suppose h{I) G Cp^. and / G Glp^o- while it is sufficient to assume 
that, for all I G Gp, fk{I,') £ C^(M'''). 

Similarly to [FW15b], we introduce the following 

Hypothesis 2.1 (Time decay). There exists Mf > 0 and a G (0,1) 

||/(/,(p,^)||^_,<M;e-“^ (3) 

Set e := iMf. We firstly state the following 

Theorem 2.2 (Strong aperiodic Birkhoff). Consider Hamiltonian (1) with h{I) := uj ■ I, under 
the Hypothesis 2.1 and the described regularity assumptions. Then, for all a G (0,1) there exists 
Ea > 0 such that the following statement holds true. For all e G (0,6^], it is possible to find 
0 < Pif < po < p and 0 < u* < uq < cr and an analytic, canonical, e—close and asymptotic to the 
identity change of variables B : 'Dp^^u^ —)• for all t G M'*', 

casting Hamiltonian (1) into the strong Birkhoff normal form 

^(oc) ^j{oo) ^ ^(oo) ^ ^(oo)) ^ ^ . j(oo) ^ ^(oo) ^ ( 4 ^ 

Hence, in the new variables, the flow with frequency w persists for all u, regardless of the 
numerical features of this vector, i.e. more specifically, no matter if it is resonant or not. The 
absence of a non-resonance hypothesis on oj implies also that (4) holds also if uj has an arbitrary 
number of zero components, implying the persistence of any lower dimensional torus. 

With a straightforward adaptation of the notational setting, the result in the general case states 
as follows: 

^In particular, if g G £p,o- then \gk\p < ||<?||p „ exp(—Ifcla) for all t G K’'". 
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Theorem 2.3. There exist e* > 0 and r G N \ {0} such that, for all £ G (0,e*] it is possible to 
find an analytic, canonical, e—close and asymptotic to the identity change of variables {I,(p,r]) = 
Afr ■ ^ ^ casting Hamiltonian (1) under the 

Hypothesis 2.1, into the normal form of order r 

if M t) = /i(lM) + (/(") ,‘P^^\t), (5) 

where 7^1'’+^) is “exponentially small” with respect to r and vanishes for'^ t +oo. Moreover, 
for all /(O) G G one has in (1): \I{t) —/(0)| < y/epo/S for all t G M'*'. 

Similarly to [FW15b] (and the mentioned previous papers), no lower bounds are imposed on 
a so that the decay can be arbitrary slow. The (natural) consequence is that either or e* 
decrease with a, see (15) and (56). 


Part I 

Proof of Theorem 2.2 


3 The normalization algorithm 


Given a function G := G{I,(p,t), define the Lie series operator exp(TG') := + 

where CgF := {T, G} = Ftp ■ Gj — Gp ■ Fj — F^Gt- The aim is to construct a generating sequence 
such that the formal limit 

B := lim B^^'> o 0...0 ^(0), (6) 

j->oo 


where B^^^ := exp(T^(j)) is such that ^8 o 77 is of the form (4). The following statement shows 
that this is possible, at least at a formal level 


Proposition 3.1. Suppose that for some j G N Hamiltonian (1) is of the form 

H^^^ = uj-1 + rj + F^^\l,ip,t). 

Then 77b+^) := B^^'l o H^^^ is still of the form (7) with 


s>l 


(s + 1)! X' 




( 7 ) 

( 8 ) 


provided that solves the homological equation 

(9) 

Since Hamiltonian (1) is of the form (7), one can set H^^^ := H with F^^'l := if. Thus, by 
induction, the form (7) holds for all j G N. Clearly, this does not guarantee that the objects 
involved in the algorithm are meaningful for all j, as it is well known their sizes can grow 
unboundedly as j increases, as a consequence of small divisors phenomena. The aim of Section 
4 (and in particular of Lemma 4.5) is to show that this is not the case: the key ingredient is the 
time decay of /. 

"‘See bound (54). 
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Proof. We get exp(£^o) = I ■uj + r} + {I, cp, t) + (w • / + r/) + Es>i (l/'SO^* o) + 

^^>2(l/s!)i2®y) (w • / + 77). The sum between the third and fourth terms of the r.h.s. of the 

latter equation vanishes due to (9). As for the last two terms, by setting as the sum of 

them, one gets = J2s->i{l/sl)C^(j)[F‘-^'> + (s + (j){uj ■ 1 + 77 )], which immediately 

yields (8) by using (9). □ 

The (formal) expansions and F^^^ = yield 

(9) in terms of Fourier components 


dtc^i\l ,t) + i\{k)c^l\l ,t) = 
with \{k) := oj ■ k. The solution of (10) is 


L Jo 


( 10 ) 


( 11 ) 


where (7,0) will be chosen later. 


4 Convergence 

The classical argument requires the construction of a sequence of nested domains C 

'Upj,aj 3 ( 7 ( 1 ), 99 ( 1 ) , 77 ( 1 )), such that Bj : Fj+i —)• Fj- The resulting progressive restriction is 
essential in order to use standard Cauchy tools, see Prop. 4.1. The estimates found in Lemma 
4.2, concerning the solution of equation (9), will be used to prove Lemma 4.5, providing in this 
way the bound on F^^'> defined in Prop. 3.1. This is achieved for a suitable sequence of domains 
prepared in Lemma 4.4 via {pj} and {ctj}. This allows us to conclude that the perturbation 
term is actually removed in the limit (6). 

The final step consists of showing that B defines an analytic map B : Fp^ 3 (7^°°^, , 77 ^°°)) —>■ 

JJpOy'^o 3 (7^ , 77 ( 0 )) = (^I^cp^r]), where p* < pj and cr* < aj for all j E N. This property is 

shown in Lemma 4.6. As will be the domain of analyticity of the transformed Hamiltonian 

via B, it will be essential to require that p*,cT* > 0. 


4.1 Some preliminary results 

Proposition 4.1. Let 7^, G : x T” x R'*' —>■ C such that and are 

bounded for some d',d'' E [0,1). Then, defining 5 := \d' — d"\ and d := max{d',7"}, for all 
d £ (0,1 — d) one has for a// s E N \ {0} 


ll'^G'^ll (1—rf—d)(p,(T) — g2 


2 e 


d{d + 5s)po- 




(l-d")(p,o-) 


\{l-d'){p,a) ’ 


( 12 ) 


where Sg = 5 if s = 1 and is zero otherwise. 


Proof. Straightforward from [Gio03, Lemmas 4.1, 4.2]. □ 

Lemma 4.2. Suppose that F^^'l satisfies exp(—at) for some M^J'l > 0, p < p 

and d < a. Define C^j := 1 + |a;|, then for all 6 E (0,1) the solution of (9) satisfies 


X 


U) 


< 

(1-5) (p,a) 


A7(J) 

a 




(i-5){p,a) 




m(^') 

a 



(13) 
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Proof. First of all, by hypothesis exp(— |/i;|(T — at), in particular, by choosing 

0) := —J^^exp{iX{k)s)fjf\l,s)ds we have that 0)| < +oo for all I E Gp- Substi¬ 
tuting c^^^(/,0) in (11) one gets \fjf\l, -3)1(13 < / a) exp{—\k\a — at) which 

yields^ the first of (13). As for the second of (13), it is sufficient to use (10), which implies, 
\dtcl^\l,t)\ < {M^^^/a){l + |a;||/c|) exp(— |/c|(7 — at) then proceed similarly. □ 

Remark 4.3. It is immediate to notice that a hypothesis of non-resonance on uj does not 
substantially improve the bounds (13). A more careful computation yields 

\c‘i\l,t)\ < + (w • A:)2)-le-l^l"^-“*, 

Hence the estimate cannot be refined due to the presence of |cq'^^(/, t)|, no matter what the 
minimum value of {u) ■ k) is. 


4.2 A suitable sequence of domains 

Lemma 4.4. Let {dj}jgN be a (real valued) sequence such that 0 < dj < 1/6. Consider, for all 
j E N, the following sequences 

Cj+i := Ka-^df^e], {pj+i,aj+i) := (1 - 3dj){pj,aj), (14) 

with K > 0 and t := 2n + 3. Then, for all 0 < pQ < p, 0 < aQ < a and cq < £a where 

£a < aA-i(27r)-2-, (15) 

it is possible to construct {djjjgpf such that (/?*,it*) = {l/2){pQ,aQ), in particular they are strictly 
positive. Furthermore limj_).oo = 0. 

Proof. Choose ej := eo(j -|- 1)“^’^ (so that hmj_^(X) = 0 by construction). By the first of (14) 
one gets 

d, = (eoAa-i)^(j + 2)V(j + l)", (16) 

hence, by (15), dj < vr“^(j -|- 1)~^. This implies — 1/^ and then, trivially, dj < 1/6 for 

all j E N. Now we have® lnnj>o(l — 3dj) = — ~ln2, hence 

\\mj^ooPj = Ponj>o(l — 3dj) > pq/2 =: p*. Analogously cr* := (Tq/2. □ 


4.3 Bounds on the formal algorithm 

Lemma 4.5. There exists K = K{pQ,crQ) > 0 such that, if e < £a where £a satisfies (15), then 


p(i) 


{pj 


< CjC 


—at 


(17) 


for all j E N. Hence, the transformed Hamiltonian B o H is in the form (4). 

^Recall (2), then use the inequality exp{—S\k\a) < Its variant X^fcg2»i(l -I- 

|ti;||fc|) exp(—dlfcla) < is used to obtain the second of (13). 

®Use the inequality ln(l — x) > —2a;In 2, valid for all x £ [0,1/2]. 
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Proof. By induction. Note that (17) is true for j = 0 setting eo := £• The condition on e ensures 
the validity of Lemma 4.4. Hence, supposing (17), by Lemma 4.2 and Lemma 4.4, we get 


X 


(J) 




\2n „ — l j—2n ^—at 


< ejiejn^Y'^a ^d- “’"e 


(18) 


By (8) and Prop. 4.1 with d' = dj, d” = 0 and d = dj (the condition dj < 1 — dj holds as 
dj < 1/6) 




{l-2dj){pj,crj) 




S>1 


< 2-^0 


(1 2dj)(pj,aj) 


pU) 




(19) 


where^ 

0 := 2ejnCui{e/a^Y< 1/2 ( 20 ) 

is a sufficient condition for the convergence of the operator exp(/l^(j)), from which X/s>i 
Hence, by (19), (20), then by (18) one gets (use also a^.,p^,,dj < 1) 




(l-2dj){pj,aj) 


<e]nC^[e/aYYp*^a ^d^e 


( 21 ) 


The latter is valid a fortiori in 'P(i- 3 dj){pj,crj)- 

In conclusion, by choosing K := nCu){e/aY)^pf^ = 2'^^^nCu,{e/a^Yp^^ from the first of (14), 
we have that (17) is satisfied for j ^ j + 1. Furthermore, by the first of (14), condition (20) 
yields 1 > AejKdja~^dJ'^ = 4dj(ej_|_i/ej)e““*. The latter is trivially true for all t G M+ by 
the monotonicity of ej and as dj <1/6. Furthermore this implies 

0 < 2dje~^\ ( 22 ) 

Hence exp{C^(j)) is well defined for all j G N. □ 

In this way the value of £« mentioned in the statement of Theorem 2.2 is determined once 
and for all. 


4.4 Estimates on the transformation of coordinates 

Lemma 4.6. The limit (6) exists, it is e—close to the identity and satisfies 

|/(oo) _ J|j^(oo) _ ^1 < (po/6)e-“^ IvpM - (^1 < (ao/6)e-“*, 

V 


(23) 


in particular it defines an analytic map B : 'Pp .,,0 
Pp,,a, for all t G M+. 


Proof. Let us start with I. Note that 


po,(To ond is an analytic function on 

< n{edjpj)-^ 


(1 2dj)(pj ,crj) 

by a Cauchy estimate (see [Gio03, Lemma 4.1]), so that the presence of n in (20) is justified. 
Hence use Prop. 4.1 with F •(— £,^(j)Yd+P^ g g — obtaining 

e“^s!0*po- This implies 




< 




S>1 


{l-3dj){pj,<7j) 


{l-3dj){pj,aj) 


< 2 ^Qpo < djPQe 


^The reason for using nCu, in the definition of 0 will be clear In the proof of Lemma 4.6. 
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by (22). In particular is e— close to the identity by (16) for all j E N, hence 

|/(oo) _ J| < - lU)\ is. It is now sufficient to recall Ylj>o^j — 1/6 in order to 

conclude. 

The argument for ip is analogous while the variable r] requires a slight modification. In particular, 
as one needs to set F •(— th® nse of the second of (13) requires the contribution 

of Cuj in (20). 

In conclusion, the obtained composition of analytic maps is uniformly convergent in any compact 
subset of Fp^ • This implies that B is analytic on Fp^ by the Weierstrafe Theorem and hence 
the image of H via B is an analytic function in the same domain. □ 


5 Further perturbation examples 

In this section we consider two alternative examples of perturbation. The main purpose is to 
show that the hypothesis of summability in time over the semi-axis is the only key requirement 
for the argument beyond the proof of Theorem 2.2. 

In particular, we shall firstly consider a decay which is assumed to be quadratic in time, while in 
the second example a perturbation exhibiting a finite number of (differentiable) bumps is exam¬ 
ined. The procedure is fully similar, with the exception of some bounds that will be explicitly 
given below. 


5.1 Quadratic decay 

Let us suppose that (3) is modified as 

In the same framework, it is immediate to show that the analogous of Lemma 4.2 yields the 
following estimates 







(1-5) (p,a) 


Clearly, in this case, the integration has led to a “loss of a power” in the decay. This is harmless 
as, by (19), a ) ~ 0{F^)0{x^^'^) + h.o.t. and then ~ (t-|-l)“^ < (t-|-l)“^ 

so that the scheme can be iterated®. 

The rest of the proof is analogous provided that the term is replaced with 1 in the remaining 
estimates. 


5.2 Differentiable bumps 

Let L E N \ {0} and h > 0. Consider an increasing sequence {ti}i=i,..,,L G such that 
— ti > 2h, then the following function 



{ai/h‘^)[{t - ti + h){t - ti - /i)]2 
0 


t £[ti- h, ti + h] 
otherwise 


®A similar (and even stronger) phenomenon could have been noticed in the original setting. Namely, suppose 
by induction that exp(—a^t). By Lemma 4.2 and (19), one finds that < 

ej+i exp{—2ajt) and so on. This leads to a remarkable rate of decay {aj = 2^a) but not to a substantial improve¬ 
ment of the estimates and of the threshold (15) of Sa, as these are uniform in j. 
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where a; E M. Considering a function E we set as 

L 

l=i 


In such case we find 


(l-5)(p,a) 



(i-S){p,a) 


with A := \^l\- The remaining part of the proof is straightforward with the obvious 

modifications. In particular, as for the proof of Lemma 4.5, one finds K = 2nC^{e/a^YhAp~^. 


Part II 

Proof of Theorem 2.3 

In order to simplify the notation, we shall use {ph,o'h) in place of {p,cr) and {p,cr) in place of 
(/5o,do) from now on. 


6 Formal algorithm 


As in [Gio03], we write Hamiltonian (1) in the form 




where 

Ho{I, p) := h{I) + p, HYI, t) := h(I, 

k€As 

where := {/c E Z”" : (s — l)iV < |A:| < sN} and iV E N \ {0} is meant to be determined. 
Given a sequence of functions —>■ C, the Lie transform operator is defined as 


Tx ■■= E 

s>0 



Id 
1 ^ 

i=i 


s = 0 
s>l 


(24) 


Let r E N \ {0} to be determined. A finite generating sequence of order r, denoted with yM, is 
such that = 0 for all s > r. Our aim is to determine it in such a way the effect of Hi ,..., Hr 
is removed, i.e. 

i/M := = Ho + n<^''+^\l,ip,t), (25) 

where the remainder contains H^r and a moltitude of terms produced during the nor¬ 

malization, which Fourier harmonics lie on A>r. The smallness of the remainder is an immediate 
consequence of the decay property of the coefficients of an analytic function. The procedure is 
standard: condition (25), with the use of (24), yields a well known diagram which s—th leveY is 
of the form 

s-l 

£s •= EgHo -|- Eg-iHi + Hg = 0, (26) 

1=1 

® Namely, those terms of the diagram which Fourier harmonics belong to A^. 
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if s = 2,... ,r and EiHq + Hi = 0 if s = 1. As sum of all the “non-normalised” levels, the 
remainder easily reads as 

7^{^+l) = ^4• (27) 

s>r 

By writing the first term of (26) in the form Eg = E^(s) + using the 

manipulation described in [GioOS, Chapter 5], one obtains a remarkable cancellation of the 
contribution of Hq. In this way, the generating sequence is determined as a solution of 




^ ;= 


Hi 


s = I 


S—1 


Hg + ^ ^ —Eg—jHj s > 2 


(28) 


A formal expansion of and of 'I's := yields for all s = 1,., 

dtc^^\l,t) + ■ k)c\^\l,t) = k G A^, 


• ,r 


(29) 


where, as usual, io{I) := dih{I). 

Remark 6.1. As a substantial difference with the isochronous case, the function oj{I) is a 
complex valued vector as I ^ Qp. In this way the exponent X{k)t appearing in formula (11) is 
no longer purely complex. More precisely, one finds a term of the form exp((a;c'(7) • k)t), having 
denoted uj{I) = U}R{I)+iuJciI), u)r^c{I) £ The size of this term cannot be controlled without 
a cut-off on k. By restricting the analysis on the levels Ag and using the fact that \ijJc{I)\ —t 0 
as p —?■ 0, a loss “of part of time decay” at each step (see Lemma 7.1) will be the key ingredient 
to overcome this difficulty. The mentioned elements are clear obstructions to the limit r —>■ oo. 


7 Convergence 

7.1 Set-up and some preliminary results 

The use of the analytic tools requires the usual construction of a sequence of nested domains. 
We shall choose, for all s = 1,..., r, the rule 

dg := d{s - l)/r, (30) 

with d G (0,1/4]. Clearly dg < d for all s = 1,..., r. Consider also the monotonically decreasing 
sequence of non-negative real numbers {««} defined as follows 

ag+i := ag{2r - s)/{2r), m := a. (31) 

Given the analyticity domain of H expressed by {ph,o'h), set a := (JhI2. Now consider the 
function n(p) := sup/gg^ |a;c'(7)|) clearly 0(0) = 0. From now on we shall suppose that p 
satisfies the following condition 

4rArO(p) < a. (32) 

The analyticity^*^ of h{I) implies the existence of Ch G [1, -|-oo) such that the value of p can be 
determined as 

p := min{pH, a{^rNCh)~^]-, (33) 

^°Obviously, Sl(p) = 0 for all p in the case of an isochronous system, so that (32) would impose no restrictions 

on p. 
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once r and N will be chosen. 

The scheme is constructed in such a way one can set (/>*, o'*) := (1 — d){p, a). 

As a consequence of Hypothesis 2.1 and of the standard properties of analytic functions, one has 


m>l, (34) 

with T := eT, where (see [Gio03, Lemma 5.2]) T := [(1 + exp(—( t/2))/(1 — exp(— it/2 ))]”' and 


h := exp(—A^( t/2). (35) 

Lemma 7.1. Suppose that ||'hs||(i_d,,)(p ct) ^ exp(—a^t), for some > 0. Then the 

solution of (28) satisfies 


4a 



(l-ds+l/2)(p>0') 




(l-rfs+l/2)(p,o-) 


(36) 


where Cr '■= /(Tf^. 

Proof. Use (29). Similarly to Lemma 4.2, we choose (1,0) := —/jg+exp(—(a;(/)-A:)r)V’fc*^ (1, r)dr. 

Note that |c^'^^(/,0)| < exp(—(1 —ds)|A:|(T) exp(|a;c'(L)||fe| — a)r)dT < +oo on by (32). 

By using again (32) one gets 

Jt o 

The first of (36) is easily recognised^^ by (31). The second of (36) follow from (37) and from 
(29). □ 

Lemma 7.2. Let A, r,r > 0 and consider the real-valued sequences {ks}s>i o-nd { 7 z}z>o defined 
as 

S — 1 I 

Ki := 7 / := T (38) 

i=i i=i 

where ki and 70 are given. Define A := r + T, then for all s >2 and I > 1 

= (Tki + rA)A^-2, 7 ; = 7orA'-^ (39) 


Proof. We shall denote with (38a) and (38b) the first and the second of (38), respectively. The 
same for (39). Let us suppose for a moment that (39a) is proven, then choose A = Tyo and 
Ki = Tyo = 7 i. By substituting in (39a) one immediately gets (39b). Hence we need only to 
prove (39a). 

For this purpose we use the well-known generating function method (see e.g. [Wil06]). Namely, 
define g{z) := Yl^=i'^nZ^, multiply each equation obtained from (38a) by as s varies, then 
“sum” all the equations. This leads to g{z) = [1 — ^z\~^{k,i{z — tz‘^) + Atz^) = (1 + Az-|-A^z^-|- 
.. .){ki{z — Tz^) -L Atz"^) = Kiz -\- (Tki -|- tA) Yln >2 which is the (39a). □ 

^^Use the inequality X^n,|>(s_i)jv — exp(—— l)cr)(^^To exp(—5mcr))" < (2/5)", where in 

this case S := — da = d/{2r). 
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7.2 Bounds on the generating function 

Proposition 7.3. For all s < r, the following estimate holds 


IIXs|l(l-d^+l/2)(p,o-) 




( 40 ) 


where the sequence {/3s}s=i,...,r G is determined by the following system 


0i 



with {0;}i=o,...,r-i G IK"*" o,nd 

r := lQnr‘^CrF{ad^pa)~^ , 
under the conditions'"^ j5i = 9 q = 1. 


(41) 


(42) 


First of all note that by (24) and (34), one has a) — ^ exp(—ait) and ||-E'o7fm||(i_rf)(p < 

jrjpn-i 0xp(—ait) (recall (31)). Hence, given by s < r, we can suppose by induction to know 
Pi,, Ps-i and 00,mj ■ ■ ■; 0s-2,m) for all m > 1 , with /3i = 1 and 0o,m = such that the the 

following bounds hold for all j = 1 , ... , s — 1 and t = 0 , ..., s — 2 

ll'^'j'll(i-di)(p,a) (43a) 

||-E'ii7m||(i_rf;^^)(p^o-) — (43b) 

By (43a) and Lemma 7.1, the bound (40) holds with j in place of s. Hence by Prop. 4.1 
with G = F = Es-j-iHm then d = maXj=i,,„,s-i{dj+i/ 2 ,4-i} = 4-1/2 and finally 
d := ds — 4-1/2 = d/ (2r), one has (by setting (5 = 0) 


Fs—j — lHm 


< Sr'^ied^pa) ^ w m ^ w \( \ 

< rj'47z_je““*+i* 

(44) 

where the property aj+i+ a;_j+i > a;+i has been used. Recalling (24), we have that (43b) holds 
also for t = s — 1, where 

r . ^ 

\m = T (45) 


I 


1=1 


Furthermore, it is easy to show from the latter that Oi^rn = in such a way, defined 

9i '■= 01,1 one gets di,m = and then the second of (41), provided 0o = 1- In conclusion, by 

using (34), and the second of (41) in the definition of iks as in (28), we get that (43a) is satisfied 
if Ps is defined as in the first of (41). Bound (40) follows from Lemma 7.1. 


Proposition 7.4. The sequence Ps defined by (fl) satisfies 


Ps < 4 V'S; 


(46) 


for s = 1,... ,r, if 


T := eh, F < h/{2r'^). 


(47) 


^^From a “computational” point of view, first compute di then proceed with Ps,ds for all s = 2,..., r. 
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Proof. The property (46) is trivially true for s = 1, hence let us suppose it for j = 1,..., s — 1 
and proceed by induction with r to be determined. Define 6i := /j, then 9i := 9i/l, 

obtaining §1 = P Clearly 9i<9i< 9i/l, furthermore 9q = 9q = 9q = 1. Hence, 

by Lemma 7.2 we have 

9i < (48) 

Now choose r,T as in (47). By using (34) and (48) in the first of (41) one gets that (46) is 
satished simply by checking that the inequality 

■= ^ + i) - 

holds true for all^^ s = 1 ,..., r. □ 


7.3 Estimates on the coordinates transformation 

From now on we shall suppose that h and e are chosen in such a way 

8 e/i < 1 (50a) 

2r^r < ^/eh (50b) 

In particular, by dehnition and by (47), this immediately implies that 

4A < 1 (51) 


As in [Gio03] it is used that, despite the generating sequence is hnite, one can use the bound 
obtained from 7.3 


X 


(«) 


< (4a)-^a.F/3,e-“’-+l^ 


(l-d)(p,(T) 

with j3s satisfying (46) for all s, as it would be, trivially, /3>r = 0. 

Proposition 7.5. Define ^)) := r [r] (/, (^, ry). Then the following estimates hold 


(52) 


I - 


(l-d)(p,o-) 


r] — r] 


(r) 


{l-d){p,a) 




(r) 


Proof. Let us start from the variable I. Firstly, note that 
In addition 


If-If 




< —e-^-+P. (53) 

(l-d)(p,cr) 8 






{l-d2)(p,a-) 


< 2nr{eda) 


-1 


X 


( 1 ) 


(l-'^3/2)(P,o-) 


{l-d){p,cr) 

< DfjTexp{—ar+it), 


(l-d)(p,a)- 


with D(j := nrCr/{2daa) by Prop. 7.3. Hence suppose a) — exp(— 0 ^+ 1 1) for 

alH = 1 ,..., s — 1 with ui = D„ and proceed by induction. 

The bound of Eil can be treated in the same way of (43b) with the difference that in this case 
the term appearing in Eil needs to be bounded separately by using (40) and a Cauchy 

estimate. This leads to ui = fiiD^ + using the same procedure used in 

the proof of Prop. 7.4 for 9i one gets Ui < /l)/2J‘~^. The required bound easily follows 

as Us < 2FDfj < Tdp < y/edp/S, where the second inequality follows from (51) and 

the last one from (50b) then from (50a). The procedure for the variables f and r/ is similar. 
The analyticity of the transformation Mr '.= T~^^ easily follows from the bounds (53) and the 
invertibility of the Lie transform operator, see [Gio03]. □ 

^^Clearly (49) holds for s < r if y{r) < exp(r — 1) for all r > 3 (let it be directly checked for r = 1, 2). Hence 
set r = n + 1 and prove that y(r)r=n+i < exp(n) for all n > 2, conclusion that is immediate as one can find that 
y{n) < n + 1 + 3e"/(4n). 
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7.4 Bound on the remainder 


Proposition 7.6. Define A := lOJ^ then for all r > 1 


■Jlir+l) 


(l-2d) (p,cr) 




( 54 ) 


Proof. Define {p',a') := (1 — d){p,a). Now recall (27) and suppose by induction, for all I = 
1,..., s — 1, m = 0,..., s — 2 with s G N 

\\^l^o\\(l-p/s)d){p',a') ^ ll-®m-f^n|| (i_(m/s)(i)(p',o-') — -^Cm,n exp( —O^+lt). 

(55) 

Indeed one can set Co,n = ^ and ei = /3i = 1 as C^(i)Hq = —4'i by (28). We stress that, 

despite based on the same computations, the argument is conceptually different from the previous 
estimates as s G (r, +oo) and the use of 6 in (12) plays here a key role. More precisely, use Prop. 
4.1 with G = F = Eg-jHo hence d" = 0 then d = d' = 6 = d{s — j)/s from which 


d = [j/s)d. This leads to 






< T{s/j)/3jes-j exp{-ar+it), implying^ that 


the first of (55) holds for I = s provided = /3s + P the latter by Lemma 

7.2. This implies |IZ)f=i- •^('® + 1)A®“^ exp(-ar+it) by using (34) and the 
trivial bound h < A. Similarly one finds ^s,n = h'^~^ A^~^, hence 


(J^e-“’'+i*)-^77('’+^) < ^(2 + s)A"-^ = 


s>r 


r + 3 


+ 


1-A 1-A2 


< 2(r+ 4)A" 


by (51). Noticing that D(i- 2 d){p,a) F the bound (54) easily follows from (51) and 

from the simple inequality (r + 4)e'’ < 5(4'"). □ 


7.5 Parameters choice and perpetual stability 


Let us discuss a possible choice of the parameters in such a way the convergence conditions are 
satisfied. More precisely by (35), condition (50a) holds if A = |'2cj“^(l + 31og2)], where [•] 
denotes the rounding to the greater integer. This implies that h > l/(16e), hence (50b) holds 
if 2^er2r < y/e. Hence, recalling (32) and (42), this condition is achieved by choosing (see also 
[GG85]) 




g2dn+2p^^2 

2‘^n+19 


(56) 


where^^ 7 = 5 + n and [-J denotes the rounding to the lower integer. The condition £ < £„, as 
in the statement of Theorem 2.3, clearly ensures that r > 1. The final value of p is determined 
with (33). 

Let us write the usual bound \I{t) — I(0)| < \I{t) — \ — I^'’^(0)| + |/^'’^(0) — /(0)|. 

The first and third term of the r.h.s. are bounded by y/edp/8 by (53). As for the second one, 
from the equations of motion , furthermore — 

eA{eda)~^ exp(—(r + a^+if)) by a Gauchy estimate and by (54). Hence 


_ jM(o)| < £A{eda)-^e-^ f e-^^+^^ds < eA{adea)-\2/ey, (57) 

Jo 

^^The use of (12) with <5 = 0 would have given [s/jfi instead of (s/j), producing in this way a troublesome 
factorial in the estimates. 

^®Note that the threshold ej takes into account of the condition (33) as we have used the obvious lower bound 
p > apH{4:rNCh)~^, immediate from (33). 
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as ttr+i = a{2r — l)(2r — 2)... (r)/(2r)^ > a2 ^. 

Remark 7.7. The bound (57) is the key element beyond the perpetual stability, despite a normal 
form of finite order. The remainder, which is bounded by a constant in the classical Nekhoroshev 
estimate and then produces a linearly growing bound for the quantity — /^'’^(0)|, is now 

summable over M'*'. Hence, a restriction to exponentially large times is no longer necessary. 

It is immediate from (57) that for all e < e* one has — /*^'’^(0)| < 2e*A{ade‘^ 

which is clearly smaller than ^/edp/A by (56). Hence \I{t) — /(0)| < y/edpj2. 
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